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We consider an effective ffamiltonian description of critical wetting transitions in systems with short- 
range forces at a corrugated (periodic) wall. We are able to recover the results obtained previously 
from a 'microscopic' density-functional approach in which the system wets in a discontinuous 
manner when the amplitude of the corrugations reaches a critical size A*. Using the functional 
renormalization group, we find that A* becomes dependent on the wetting parameter ua in such a 
way as to decrease the extent of the first-order regime. Nevertheless, we still expect wetting in the 
Ising model to proceed in a discontinuous manner for small deviations of the wall from the plane. 



I. INTRODUCTION 

In a recent article ^ we have studied wetting transi- 
tions at a periodic non-planar substrate (i.e. a corrugated 
wall) within a Landau-like mean-field (MF) square gradi- 
ent theory, which is appropriate for modelhng adsorption 
in fluid and simple magnetic systems with short-range 
forces . The surprising conclusion of this analysis was 
that second-order (critical) wetting transitions occurring 
in the planar system are generically corrugation-induced 
first-order when the root mean square amplitude A of 
the corrugations exceeds a tricritical value A*, whose 
numerical value is less than a bulk correlation length. 
The analysis of the Landau model is based on a 
perturbative expansion of the free energy about the 
equilibrium value of the planar system. It transpires 
that the free energy correction due to corrugation can 
be expressed in terms of (planar) correlation functions, 
and is amenable to a graphical interpretation which 
generalizes the well known Cahn construction |^ for the 
planar system. Despite the convenience of a geometrical 
description and the simplicity of the final results, the 
calculation is rather involved and perhaps obscures the 
underlying mechanism for the shift in the order of the 
transition. In this paper, we seek to elucidate the 
simplest possible effective Hamiltonian theory consistent 
with the previously derived results for the shifted MF 
phase boundary and proceed to discuss the influence 
of thermal fluctuation effects beyond MF level on these 
predictions. As we shall see this is a surprisingly subtle 
problem and complete answers cannot be given for the 
renormalized phase boundary. Indeed, even at MF 
level problems arise since a simple effective Hamiltonian 
approach does not quite recover all the quantitative 
results of the Landau theory. In discussing these 
issues we also seek to explain why a previous effective 
Hamiltonian study B did not find any evidence for a 



corrugation-induced first-order wetting transition in the 
non-planar geometry. 

In the last ten years or so one of the main thrusts 
in theoretical wetting research has been the study 
of adsorption in such non-planar geometries. The 
conformation of thin liquid films on rough surfaces has 
been described quite extensively Q while the influence of 
disordered (self-affine) substrates on three dimensional 
wetting transitions has also been investigated, using 
both replica and renormalization group methods 
New critical behaviour is indicated when the roughness 
of the substrate exceeds the thermal roughness (as 
measured by the roughness exponent 0|) of the wetting 
layer. However, the transition remains second-order in 
character. In two dimensions, this no longer need be true 
and even richer behaviour occurs Work has also be 
carried out for systems with van der Waal's interactions 
for both the wedge ^ and groove geometries and 
quite general conditions for roughness-induced wetting 
have been found . 

To begin we recall our basic 'microscopic' model of 
fiuid adsorption at a corrugated wall (appropriate to 
systems with short-range forces) and discuss the results 
and interpretation of our earlier work. This is of some 
importance, as we shall show that the simplest available 
phenomenological approach does not quite agree with all 
our previous predictions. 



II. PERTURBATIVE ANALYSIS OF A LANDAU 
THEORY 

Writing zw{y), with y = (^1,2/2), for the local 
height of the wall above the z = plane we consider 
a Landau- Ginzburg- Wilson (LGW) Hamiltonian of the 
order parameter (magnetization) m{v), 
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HLGw[m] ^ dy< dz 



Mmi) \ (2.1) 



where mi denotes the surface magnetization at vector 
position y along the z — plane, i.e. mi{y) — 



m{y,zw{y)), while c/){m) and <j)i{m) 



him 



are appropriate bulk and surface potentials. Here c and 
hi are the surface enhancement and field, respectively. 
The final multiplicative term in (2.1) represents the local 
increase in area due to the corrugations of the wall 
expanded to second-order in zw- Note, we have written 
V = (Vj_, -§;) with = We suppose that 

the bulk exhibits two-phase coexistence in zero bulk field 
(h = 0) between phases with magnetization nia (> 0) 
and (< 0) for T < Tc- Furthermore, we set ft-i < 
and focus on wetting of the wall-a interface by the (3 
phase. 

At mean-field level, corresponding to Landau or square 
gradient theory, the Hamiltonian needs to be minimized 
to obtain the equilibrium magnetization. This is a 
straightforward task for the planar system zw = 0. 
Denoting n = (f>"{mfj), the inverse bulk correlation 
length of the /3 phase, we recall that for c greater 
(less) than k the wetting transition is second (first)-order 
Q. For the non-planar system, however, even the MF 
calculation is non-trivial due to the loss of translational 
invariance along the wall. Previously [0, we have 
developed a perturbative approach to the problem in 
which the equilibrium free energy of the wall-a interface 
is written 
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2(27r)(''-i) 



dq g^A^(q)|zH'(q)|' 



(2.2) 



Here the first two contributions represent the bulk 
and surface free energy (tension) of the appropriate 
planar system while the final term is the non-planar 
correction written perturbatively in terms of the Fourier 
amplitudes of the wall function zw{y)- The quantity 
to be determined A^(q) has the dimensions of a surface 
tension and by construction satisfies the long wavelength 
identity 



A7r(0) = (7wa 



(2.3) 



required from infinitesimal rotational invariance. More 
generally, A^(q) can be related to the surface magneti- 
zation TOtti and correlation functions of the planar system 



9^A^(q) = (^i(m^i) 
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G(0,0;q) G(0,0;0) 



(2.4) 



Here (5(^1, 2:2; q) = G{zi,Z2;q) is the transverse Fourier 
transform of the connected two-point correlation func- 
tion, G(ri,r2) = (m(ri)m(r2)) - (m(ri))(m(r2)), for 
positions distances zi and Z2 from the wall, i.e. it is 
surface correlations which determine the free energy 
correction. In this way one may derive an elegant relation 
for A^(q) (= A^((7) from (2.4)) which complements the 



well known expression for the surface tension 



(Jwa = (piim^i) + / dmQo{m) 



(2.5) 



where = TO7r(oo) is the bulk magnetization and 
Qo{m) — ■\/2(0(m) — (f)(mb)) is the usual function 
appearing in the Cahn construction for planar square 
gradient theories. The final expression for A^(q) is 

A^(9) = 0i("^7ri) + / dmQ{m;q) (2.6) 

where Q{m; q) is a modified Cahn function satisfying the 
differential equation 



d 
dm 



= q^Q 



(2.7) 
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Clearly, when q = we have Q{m; 0) 
reduces to awa as quoted earlier. 

Application of this perturbation theory to the case of 
adsorption at a corrugated sinusoidal wall. 



zw{y) = V2Asm{pyi) 



(2.8) 



with A the root mean square width of the corrugations 
and 27r/p their period, yields results best expressed in 
terms of the contact angle 6^^ = 0-^{T, c, hi) of the planar 
system. Thus, for the case of strongly first-order wetting 
transitions in the planar system occurring at temperature 
T^, say (for which 9^ (x {T^-T)^), it was found that the 
corrugated wall-a interface was completely wet by the (3 
phase at a lower temperature satisfying 



pA 



for c < K 



(2.9) 



(valid for p <^ k) with the transition remaining first- 
order. Indeed, the possibility of corrugation-induced 
second-order wettin g tra nsition s is ruled out completely 
from the expansion (2^) with (2.4). 

On the other hand, for planar second-order wetting 
transitions (for which 0^ oc (T^ ~ T)) the corrugated 
geometry showed a first-order phase transition for suffi- 
ciently large ^ > A* at a reduced temperature satisfying 



e,^p[A^-A*\^ 



for c > K 



(2.10) 



For A < A* the transition remained second-order and 
occurred at the same planar wetting temperature T^.. 
Surprisingly, the threshold or tricritical amplitude A* is 
comparatively small and is explicitly determined as 
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kA* 



for c> K 



(2.11) 



which is less than a bulk correlation length. Thus, 
even minor deviations from the plane can lead to a 
corrugation-induced first-order wetting transition. 

Before trying to rederive these results using an effective 
Hamiltonian approach it is worth mentioning a few points 
concerning their interpretation. Firstly, the first-order 
result (2.9) is precisely the expression obtained from 
a naive application of Wenzel's empirical law ||l^ to 
wetting transitions. Recall that Wenzel observed that 
the contact angle 9p of a droplet on a rough surface (of 
area Ap) appeared to satisfy the relation 
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A. 
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(2.12) 



For a corrugated wall 
order in pA) and setting 
small Ott. However, the macroscopic Wenzel law (|2.12 



^1 (to q uad ratic 
recovers ( 2S ) fo r 



is certainly not universally valid and is restricted to 
adsorption problems in which the transverse correlation 
length characterizing the capillary wave-like fluctuations 
of the a/3 interface is not much larger than the bulk 
correlation length. This condition is met for strongly 
first-order phase transitions since the thickness of the 
wetting P layer does not exceed a few bulk correlation 
lengths before the transition to infinite adsorption occurs. 

Secondly, for c > k (i.e. planar second-order wetting 
transitions) the amplitude A* vanishes smoothly as the 
planar tricritical point c = k is approached. We 
emphasize here that this is required in order that the 
global surface phase diagram shows a smooth cross- 
over to the Wenzel-like result ( |2.9| ) appropriate for 
the planar first-order regime. Given that we can 
discount the possibility of corrugation-induced second- 
order behaviour in the c < k sector, a non- vanishing value 
of A* as c — > would somewhat surprisingly result in a 
discontinuo us su rface phase diagram. Consequently, the 
prediction ( ^.11 ) imbues the surface phase diagram with 
a natural topology facilitating a smooth cross-over near 
planar tricriticality. Sections through this diagram are 
sketched in Fig. |l| and Fig. ||. 

Finally, we note that deep in the planar second-order 
sector c ^ K the result for A* exhibits universal-like 
properties 



kA* 



1 



for c ^ K 



(2.13) 



This prediction is appropriate for planar critical wetting 
transitions close to the bulk critical temperature and 
sufficiently far from the planar tricritical point. As we 
shall see, the right hand side of this equality is associated 
with the numerical value of a hyperscaling amplitude. 
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FIG. 1. Schematic wetting phase diagram for fluid adsorp- 
tion in a system with a corrugated wall and c > k. The 
lines Cl and C2 are loci of first and second-order wetting 
phase transitions, respectively, which meet at the tricritical 
point corresponding to A = A* . The vertical axis is a linear 
measure of the temperature scale Ttt — T for small 9,^ , where 
Ttt is the wetting temperature in the planar system. 
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FIG. 2. For a system with c < k only first-order transitions 
take place. The wetting temperature is reduced below T,r by 
an amount proportional to the wall amplitude squared (at 
least for small A). 



III. EFFECTIVE HAMILTONIAN THEORY 

A. Predictions of a simple model 

The simplest continuum interfacial Hamiltonian model 
of fiuid adsorption at a non-planar wall is 

H[i,zw]~ J dyi^^j:^p{V^ef + W{e-zw)^ (3.1) 

where E^^ is the stiffness coefficient of the unbound a(3 
interface, W{£) is the binding potential for adsorption 
at a planar wall and £{y) is a suitable measure of the 
local thickness of the wetting film. This is certainly 
a plausible starting point for investigations and has 
been employed by numerous authors ||l^] to understand 
interfacial fluctuation effects at rough and self-affine 
walls. Therefore, it is surprising to note that the model 
does not recover all the Landau MF results quoted 
earlier. Nevertheless, it does manage to capture the 
correct physics away from the planar tricritical point and 
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provides a simple means of understanding the origin of 
the corrugation-induced first-order behaviour for c 3> 
K. Before discussing the effect of wall corrugation on 
planar first- and second-order wetting transitions we 
make a few general remarks. Firstly, at MF level 
the problem reduces to finding the global minimum of 
H[£,zw]- A wetting phase transition occurs when the 
global minimum approaches zero from below so that the 
free energy is identical to that of an unbound interface 
H[oo, zw] = which always represents a minimum. This 
comparison of free energies appears not to have been 
adopted by Rejmer and Napiorkowski Q (RN) in their 
effective Hamiltonian study and as a consequence they 
concentrate solely on second-order behaviour (note these 
authors actually consider a more complicated effective 
Hamiltonian model systematically derived from the LGW 
functional). Secondly, following standard theory, the 
MF equilibrium values of the planar film thickness It^, 
correlation length ^j| and singular free energy f^^J^ are 
identified through 

W'ii^) = ; C|? - TTO ; /ir„l = (3.2) 



The latter quantity is conveniently related to the planar 
contact angle by Young's equation ||^ 



fshl = <Jaf){cOs9^ - 1) 



(3.3) 



where cTc^ is the surface tension of the a/3 interface. 
Denoting £p{y) the collective coordinate field that 
minimizes H[£,Z]y] we have 

^o^p^lSip = W'{6£p ~ zw) (3.4) 

which has to be solved perturbatively. Assuming that 
kw(y)| and 5£p = £p — £tt are small, we find to first order 

S„/3Vi«p - W"{£^) (5ip - zw) (3.5) 

which is trivially solved on introducing the Fourier 
transforms 



S£p{q) = J dyc^^-yd£p{y) 
zw{<l) = J dyc'^i-yzwiy) 



We find 



.(q) = 



,2^2 
1 Sll 



1 



(3.6) 
(3.7) 

(3.8) 



and see that the 'healing length', in the terminology of 
is equal to the transverse correlation length of the 
planar system. Short wavelength surface undulations, 
such that q£_ii » 1, are completely damped by the 



interfacial stiffness. The interface is flat and the effects 
of surface roughness are negligible. We also point out 
that it is possible to go beyond the simple form of the 
effective Hamiltonian given in ( |3.l| ) and to include a non- 
local interaction between the wall and the a/3 interface. 
Equation ( |3.8| ) is then modified by the appearance of a 
kernel function which is simply unity for our case (see |^ 
for more details). 

Thus to quadratic order in zw{(l) (equivalent to the 
free energy expansion (2.2) of the Landau theory), the 
singular contribution to the surface free energy for a 
bound (non-wet) phase at a corrugated wall is 



-^non— wet — -^[^p: ^w] 

= A^W{£^) 
1 



2(27r)2 



dq 



1 + q^^n 



l^w(q)|' 



(3.9) 



which must be compared with 



Fwet = H[oO,Zw] 

= 



(3.10) 



If we specialize to the case of a corrugated wall ([2.^), 
the integral can be evaluated and the wetting phase 
boundary is then given by the solution of 



\wie„)\ 



2(i+p2ep 



(3.11) 



1. First-order planar wetting transitions 

At a first-order transition the singular free energy /sing 
vanishes linearly as T ^ and the correlation length 
^11 remains finite at the transition. Thus, provided the 
wavelength of the corrugations is much larger than 
the transverse correlation length of the (planar) thin film, 
we find a shifted first-order wetting transition occurring 
at a lower temperature 



Ap 



forp^ll<l (3.12) 



identical to (|2.9| ) for isotropic fiuid interfaces (and recall 
~ (T^ — T)^). Close to the planar tricritical point 
c = K, where ^|| is large for the th in fil m phase, this 
prediction is no longer accurate and ( 3.11 ) suggests that 
the expression for the shifted phase boundary shows 
cross-over to cx ^ for fixed p. However, we will not 
dwell on this since, as we shall show below, the simple 
interfacial model is somewhat unreliable close to planar 
tricriticality (c « k) . 
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2. Second-order planar wetting transitions 



The simple Wenzcl-likc result ( 3.12| ) is clearly inappli- 
cable for describing the effect of corrugation on a planar 
second-order wetting transition due to the divergence of 
the correlation length as T ^ T~. Instead, directly 
from (3.11), we have 



2„2 



(3.13) 



for the shifted phase boundary corresponding to a 
corrugation-induced first-order transition. In fact, the 
absence of real solutions (for 9^) in (3.13) indicates that 
the transition remains second-order. This expression 
further simplifies because the hyperscaling amplitude 
ratio 



7^ = T^T„ 



■> sing> I 
(Sa/j/K^) 



(3.14) 



is a pure number in MF theory (see below). Accepting 
this for the moment, we conclude that the interfacial 
model predicts a shifted wetting phase boundary of the 
form 



O-a/3 ^ 



for A> A* 
for A<A* 



(3.15) 



and recall ^ (T^ — T). Moreover, the tricritical 
amplitude A* is simply determined by the value of the 
amplitude ratio TZ, 



kA* = Vm 



(3.16) 



At this stage the results of the interfacial model appear 
very similar to those of the Landau theory calculation 
for which we can identify = T^ap. However, 

the prediction for A* does not quite agree as seen by 
calculating the value of TZ. If we adopt the standard 
form for the binding potential in zero field {h = 0) 



w{i) 



— ae 
oo 



-Ki 



for £ > 
for ^ < 



(3.17) 



with 



TZ — \ so that 



T-^ — T and 6 > 0, a simple calculation yields 



kA* 



1 



(3.18) 



instead of the result ( p. 11 ). Unlike the direct analysis 
of the Landau theory functional, the simple interfacial 
model predicts that the amplitude A* is independent 
of the surface enhancement c and does not vanish as 
the planar tricritical point (c = k) is approached. 
Nevertheless, this approach does agree with the results 



of our earlier theory for strongly first- and second-order 
(planar) wetting transitions corresponding to c <C k 
and c 3> K, respectively. In particular, the 'universal' 



result (2.12) can be traced to the numerical value of the 
appropriate ratio TZ of hyperscaling amplitudes. 

Before discussing the possible reasons for the dis- 
crepancy between the Landau and interfacial model 
calculations close to planar tricriticality we focus on 
the strongly second-order sector c ^ k and enquire 



how results (3.15) and (3.16) are modified by thermal 
fluctuations in three dimensions. 



B. Fluctuation effects away from planar tricriticality 

Wall corrugation has only a minor influence on wetting 
transitions below the upper critical dimension where 
entropic fluctuation effects dominate the unbinding 
mechanism |l^ . However, at the upper critical dimension 
d = 3 (restricting our attention to systems with short- 
range forces) we can anticipate that many of the qualita- 
tive features seen in the MF calculation retain relevance 
even after thermal fluctuations are allowed for. If we 



simply assume that the phenomenological model (3.1) 
is a reasonable description of fluctuation effects {away 
from the planar tricritical point) it is straightforward to 
develop a linear functional renormalization group (RG) 
analysis along the lines formulated by Fisher and Huse 
(FH) ||l^ for the planar problem zw — 0. In fact, 
assuming that the renormalized Hamiltonian H^^^ [£, zyy] 
is of the same functional form as (3.1), then the RG 
transformations are unchanged and 



zw\ 



dy { -I]„;3(Vi£)2 + e-2*Ty(*)(^ - zw) 



(3.19) 



in d = 3, where 6 = e* is the usual spatial rescaling factor 
and 



w'^*\e-zw) 



J2t 



di'W{i')(i 



Here uj is the wetting parameter 

UJ = 



(3.20) 



(3.21) 



with ks the Boltzmann constant. Following FH we 
adopt a matching procedure and choose t ^ t* such 
that the renormalized curvature (£) is about unity 

(see appendix) at its minimum. The structure of 
the perturbative analysis of (3.19) is then identical to 
that described earlier for MF theory (equations (3.4)- 
( |3.1l| )). Consequently, the expression for the shifted 
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phase boundary and threshold corrugation amphtude are 
identical to those quoted earlier (see ( ^.15 ) and ( 3.16 ), 
respectively) although the temperature dependence of 
and the numerical value of TZ are now dependent on uj. 
Thus, the planar contact angle vanishes like 



(T, - Tpi 



where [[L5| 



oo 



for <uj < ^ 
for i < < 2 
for cj > 2 



(3.22) 



(3.23) 



is the standard correlation length critical exponent 
describing the divergence of in the planar system as 
T — !■ . Furthermore, one can determine the value 
of the amplitude ratio 7?,'*-' (see appendix) so that the 
renormalized expression for the threshold amplitude A* 
satisfies 



kA* 



1 

(2a;) ^ 



for < < i 
for w > i 



(3.24) 



which we emphasize must be regarded as a 'high' 
temperature prediction valid away from the planar 
tricritical point. Thus, the effect of thermal fluctuations 
is to increase the threshold amplitude A* so that it is 
more difficult to produce a corrugation-induced first- 
order transition. This is supported by the exact results 
found in d = 2 where the transition is always second 
order independent of the size of A Nevertheless, 
the numerical value of A* remains relatively small for 
physically relevant estimates of w. For example, close to 
the bulk critical temperature u approaches a universal 
value ujc ~ 0.77 Q which implies that kA* w 1.1 if 
is close to Tc- This is hardly different from the MF 
prediction and could be checked in Ising model simulation 
studies. 

Finally, in this subsection we note that in deriving 
these results we have taken as our starting point the 
simplest phenomcnological model of fluctuation effects 
at wetting transitions and have ignored the possibility 
of including a position dependent stiffness coefficient fl^ ] 
and also coupling to order parameter fluctuations near 
the wall ||l8|-pl[|. While both these amendments to the 
standard model have some important consequences in 
d = 3 they do not effect the location of the planar 
wetting phase boundary to any great extent. Thus, 
while the current theoretical expectation is that the 
planar wetting transition (occurring in the Ising model 
say) is extremely weakly first-order beyond MF level 
(even for c > k) there is no direct evidence for this in 
Ising model simulation studies which certainly appear 
to show second-order behaviour p3[ |. Consequently, we 
feel justified in ignoring the possibility of a fluctuation- 
induced first-order transition since it is unlikely to 
interfere with the mechanism for corrugation-induced 
first-order wetting which is our central concern here. 



C. Remarks on improved effective Hamiltonians 

To complete our article we return to the discrepancy 
between the perturbative Landau theory and simple 
interfacial model predictions for the tricritical amplitude 
A* (equations (2.11) and ( 3.1^ ), respectively). Both 
results agree for c ^ k and are consistent with 
independent numerical minimization of the Landau free 
energy functional in this limit |2^ . However, they are 
qualitatively different for c ~ k since only within the 
perturbative Landau theory calculation does A* vanish, 
as it must, as planar tricriticality is approached. For, 
as mentioned earlier, a non- vanishing A* would imply 
a discontinuous cross-over to the Wenzel-like behaviour 
expected in the first-order sector c < k. 

In principle, it is of course possible to recover all the 
Landau theory using an effective Hamiltonian. If we 
denote ^(y) as the surface of fixed magnetization = 
then a constrained minimization of the LGW functional 
subject to the crossing constraint m(y,^(y)) = (see 
]l7[) defines a Hamiltonian 



H[£,zw] = mm HLGwim] 



(3.25) 



which identically must recover the MF free energy (2.2) 
on further minimization with respect to the collective 
coordinate £ {y). However, the constrained minimization 
involved in (|3.25 ) is, o f co urse, at least as difficult as the 
original MF theory of (2.1) unless various approximations 
are used [ p^ . It is these very approximations which lead 
to the discrepancy in A* described above. 

In fact, one should not be surprised at the limited 
domain of validity of ( ^3.l| ). As is now appreciated even 
for a planar substrate, the simple interfacial model does 
not provide an accurate description of magnetization 
correlations near the wall |2^,|2^ but it is these very 
correlations that are explicitly incorporated into the free 
energy correction kernel A-rr{q). To derive the correct 
Landau expression for the wall correlation function 
G(0, 0; q) a two-field Hamiltonian H[X, i] ||2l[] is required, 
where X{y) is the collective coordinate most suitable for 
modelling magnetization fluctuations near the surface. 
Note, that this need not be interfacial-like and indeed is 



not so at the critical wetting transition 19 2l[|. 

This suggests that the effective Hamiltonian for wet- 
ting at a non-planar wall should resemble the two-fleld 
models of wetting at a planar surface, provided X{'y) 
is chosen to be an interfacial-like variable, describing 
translations of a contour of flxed magnetization close to 
the value at the wall. This is indeed the case and the form 
for H[£, zw] derived by RN is almost identical to the two- 
field Hamiltonian of Boulter and Parry |l^. Using the 
notation of the latter the generic form for zw] in the 
long wavelength limit is 



H\l, Zw] = 



6 



+ ASi2(f - Zw)V_L^.V_LZW 

+ ^^22{e - zw)iy±if + W{i - zw) \ (3.26) 



E^jy constitute the elements of a stiffness- 



where the 

matrix describing the position dependent corrections to 
the separate surface tensions 



Sll = Uwa + AI]ii(i? — Zw) 
S22 = + AE22(^ - Zw) 



(3.27) 
(3.28) 



All the AYj^^{t) vanis h as £ cxd and if this position 
dependence is ignored ( 3.26| ) becomes simply (^]^) . From 
the explicit expressions for the various functions, RN 
derive the following asymptotic expansions (in zero field 
h = 0) 

AEn(£) uioX + (W20 + U2ik£)^' + • • • (3.29) 
AEi2(£) = wwX + (W20 + W2inl)X^ + ■■■ (3.30) 
AE22(£) - (r/10 + ryiiK^)^ + (7720 + Tl2i^^l)X^ + ■■■ (3.31) 
T/F(^) = vwX + w2oX2 + . . . (3.32) 

where X = e^''^ and the ellipses denote terms of cubic 
and higher order in X. 

Here we point out, for the first time, that the four 
functions /^Yj^^{t) and W{i) are not independent but 
necessarily satisfy a functional stiffness-matrix binding 
potential relation 

ASn(£) + 2AEi2(£) + AS22(£) = W{1) ~ IW'{1) 

(3.33) 



which is the analogue of the identity ( |2.3| ) in the Landau 
perturbation theory. This relation is valid for arbitrary 
choices of the potentials (j)[m) and 0i(m), and if zw] 
is defined via a partial trace [0 rather than just a MF 
saddle point identification (as in (3.25) and Q). The 
identity follows from simply requiring the Hamiltonian 
( |3.26 ) to be invariant with respect to infinitesimal 
rotations of the plane of the wall. For the equilibrium 
planar MF position satisfying W(£^) = 0, the 
identity reduces to the (bare) stiffness-matrix free energy 
relation pertin ent t o the two-field theory of wetting 
po| . Equation ( 3.33 ) is a more general requirement and 
significantly constrains the beh aviour of the coefficients 
in the asymptotic expansions ( 3.29 - 3.32 ). Specifically, 
we obtain the following 



vw = 27711 

f 10 = wio + uio + 27710 
V20 = W20 + W20 + 27720 
2w20 = W21 + U21 + 27721 



(3.34) 
(3.35) 
(3.36) 
(3.37) 



All these relations are obeyed by the coefficients explicitly 
calculated by RN using the double parabola approxima- 
tion pa], but are equally valid for arbitrary 0(m). 



Recall that a two-field approach seems sensible due to 
the appearance of the planar wall correlation function in 
( ^.4[ ). As pointed out in this correlation function 
can be recovered (in MF theory) by using a H[ii,i\ 
Hamiltonian, with two interfacial fields modelling fluc- 
tuations at the wall and at the a/J interface, respectively. 
However, such a tack is best suited for the complete 
wetting regime and at critical wetting £1 should no 
longer have any interfacial-like component p9 21| and 



as such a H[X,£] Hamiltonian is the optimal choice. 
Unfortunately, the technique of 'freezing' the lower field 
into the configuration of a rough wall (as used in ( p. 26 )) 
in order to describe wetting in a non-planar system 
is no longer valid. The field X{y) is not interfacial- 
like and consequently does not have the dimensions of 
length. Progress can be made by introducing two fields 
at the wall, one interfacial-like so that it can take up 
the wall configuration and the other non-interfacial-like 
to allow the Hamiltonian to recover critical wetting wall 
correlation functions in the limit of zw going to zero. 
However, pushing this three-field Hamiltonian beyond 
MF theory seems, at the moment, to be prohibitively 
difficult. 



IV. DISCUSSION 

In this paper we have rederived recent predictions for 
corrugation-induced first-order wetting transitions using 
an effective Hamiltonian (with a harmonic approximation 
to the binding potential) . We have shown that the phase 
diagram (Fig. |^) is valid in three dimensions and that 
the critical amplitude A*, for which walls corrugated 
with A > A* discontinuously wet, depends on the 
wetting parameter uj. The value of A* is related to a 
ratio of hyperscaling amplitudes — a prediction which is 
open to investigation by Ising model simulation studies. 
Consistent with exact results in d = 2 jlj], it can be 
seen that fluctuations extend the size of the second-order 
regime beyond the MF predictions. This extension is not 
very significant for the three dimensional Ising model at 
temperatures close to Tc- 

Importantly, our results are not quantitatively reliable 
near planar tricriticality and we are only confident of the 
predicted values of A* for wetting transitions away from 
the region c ~ k. 

We discuss the form of an improved effective Hamilto- 
nian as considered in |^]. We point out (as do RN) the 
similarity with the two-field theory of coupling effects at 
planar wetting transitions and indicate that an accurate 
description of surface correlations is needed for a global 
prediction of A*. Our analysis is consistent with that 
of RN who did not compare the bound state free energy 
with that at ^ = 00 and so did not explore the first-order 
regime. Nevertheless, they did show that any continuous 
divergence of the interface thickness necessarily occurs at 
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the same wetting temperature as in the planar geometry. 
This is entirely in keeping with our earlier remark that 
for A < A* the transition remains second-order. Using 
a new stiffness-matrix binding potential relation we are 
also able to prove that the coefficients appearing in the 
Hamiltonian, via the asymptotic expansions of AE^^ 
and W(£), are not independent but obey simple linear 
relations. 

To conclude, we believe that corrugation-induced first- 
order wetting is present in three dimensional systems 
with short-range forces but that more work is required 
to elucidate the behaviour of the threshold amplitude A* 
near the planar tricritical point. 



APPENDIX: CALCULATION OF THE 
RENORMALIZED AMPLITUDE 7^ 

As there is no non-trivial fixed point of the renormal- 
ization group representing the wetting transition, FH 
adopt (as mentioned previously) a matching procedure 
to determine the value of t up to which renormalization 
takes place. This value t = t* is chosen to be 
that at which the transverse correlation length of the 
renormalized binding potential is of the same order as the 
non-critical bulk correlation length k~^. One requires 



de 

a2^(t*)(£) 



dp 







(Al) 
(A2) 



ay 2b 2 



(A6) 



implying, via (3.16), that A* not shifted from its mean- 
field value. 



Regime II ^ < lo < 2 

The algebra is now a little more involved, 
renormalized potential is approximately 



-2t 



t^(*)(£) 



K 



e 



with K a constant |15|. Imposing ( |AlD we find 

IK 



ae 



ujt-l 



while 



e-2*W^(*)"(£)«-ae'^*-^ 



Consequently, 



making use of (A8). However, FH show that 



The 



(A7) 



(AS) 



(A9) 



(AlO) 



which defines the equilibrium wetting layer thickness i 
and allows a mean-field type an al ysis t o be appro priate. 
By setting k = 1 , equations ( |3.2|) , ( |3.14 ) and ( 3.19 ) imply 
that 



V(*)"(£) 



(A3) 



and so, similar to |15[ , we have to consider three separate 
regimes for the three different possible forms of the 
binding potential. 



and so 



V8w \t~ - \nt 



UJ 1 

21 

OJ 

2 



Int 
8t 



for large t, i.e. T close to T^. 



(All) 



(A12) 



Regime I < i 



For this case FH find that 

e-2*w^(*)(^) = -ae"*-^ + 6e4'^*-2^ 

and so 

-3ujt 



a 



2b 



From (A2), calculating TZ^^^ is straight-forward 



(A4) 



(A5) 



Regime III uj > 2 

For this case we adopt the notation of and write 
£ ~ fit, the binding potential has the form 



(A13) 



where 



-a 



b 2ujc „ , , , 
- + :T3ir + — + 0(1A) (A14) 



1 - -i^- 2 
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The constant c is introduced as the hnear renormaHzation 
group cannot handle correctly a completely hard wall, 
see ( |3.17 ), and so a 'soft' approximation is used, with 
W{£) = c for £ < 0. Note that this problem does not 
arise in certain non-linear formulations [Q. 
The matching procedure leads to ]l5[| 



£ w VSujt 
8uj) = 



(A15) 
(A16) 



where the latter defines the renormalized wetting tem- 
perature. To leading order 



-— 1 



/i X dKt 



implying that 



Writing r 



a - ac 



FH find 



0{l/t) 



(A17) 



(A18) 



(A19) 



Equation (A_f4) can then be used to show that for ^k, jl 

(A20) 



Kt = 0{t) + 0{l/t) 

K[ = constant + 0(t) + 0{l/t) 



Hence, if VS^^ is written as 



(t)2 _ 



2ujt 



(A21) 



then the term in square brackets must be almost 
constant, let it be C, s ay, for T near T^, that is r <C 1 
and t » 1. Using (|A15|) , 



7^(*) « c 



(A22) 



with C an unknown constant. However, by continuity of 



T?.*^*^ at w = 2 and comparing with (A12) we can see that 
C= 1. 
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